Measurement of Nuclear Transparency for the A{e, e'n^) Reaction 
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We have measured the nuclear transparency of the A{e,e''K^) process in ^H.i^C, "Al, "^^Cu and 
^^^Au targets. These measurements were performed at the Jefferson Laboratory over a four momen- 
tum transfer squared range = 1.1 to 4.7 (GeV/c)^. The nuclear transparency was extracted as 
the super-ratio of {cta/o-h) from data to a model of pion-electroproduction from nuclei without n-N 
final state interactions. The and atomic number dependence of the nuclear transparency both 
show deviations from traditional nuclear physics expectations, and are consistent with calculations 
that include the quantum chromodynamical phenomenon of color transparency. 

PACS numbers: 25.30.Rw, 24.85. +p 



In the context of perturbative Quantum Chromo Dy- 
namics (QCD), Brodsky and Mueller [l[ predicted that 
at sufficiently high momentum transfers, the quark-gluon 
wave packets of hadrons can be produced as a 'color neu- 
tral' object of a reduced transverse size. If this compact 
size is maintained for some distance in traversing the nu- 
clear medium, it would pass undisturbed through the nu- 
clear medium. This is the so-called phenomenon of color 
transparency (CT). Nuclear transparency, defined as the 
ratio of the cross section per nucleon for a process on a 



bound nucleon in the nucleus to that from a free nucleon, 
is the observable used to search for CT. A clear signature 
for the onset of CT would involve a rise in the nuclear 
transparency as a function of Q^. Later works f3| have 
indicated that this phenomenon also occurs in a wide va- 
riety of models which feature non-perturbative reaction 
mechanisms. 

More recently, CT has been discussed in the context of 
a QCD factorization theorem derived for meson electro- 
production which states that the meson production 
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amplitude can be expressed in terms of a hard scattering 
process, a distribution amplitude for the final state meson 
and a parametrization of the non-perturbative physics 
inside the nucleon known as Generalized Parton Distri- 
butions (GPDs) [4]. Factorization is expected to be valid 
for > 10 (GeV/c)^, however, under certain conditions 
it may also be applicable at lower It has been 

suggested that the onset of CT is a necessary (but not 
sufficient) condition for factorization to occur 6]. The 
underlying assumption is that in exclusive "quasielastic" 
production, the hadron is produced at small inter-quark 
distances. Thus, the unambiguous observation of the on- 
set of CT is critical as a precondition to the validity of 
factorization in meson production, and because it would 
open a new window to study the strong interaction in 
nuclei. 

A number of searches for CT have been carried out in 
experiments using the 2p), A(e,e'p) reactions, co- 
herent and incoherent meson production from nuclei and 
pion photoproduction reactions - [13. At high en- 
ergies the di-jet experiment at Fermi Lab and 
production at DESY Jllj are consistent with CT, and 
it is necessary to include CT to understand shadowing 
in nuclear deep inelastic scattering However, the 

predicted onset of CT at intermediate energies, typical 
of experiments at JLab and BNL, has not yet been un- 
ambiguously observed. The most recent experiment at 
JLab to look for CT in qqq hadrons using the A(e, e'p) 
reaction does not show any increase of the nuclear 
transparency up to Q"^ — 8.1 (GeV/c)^ and rules out sev- 
eral models predicting an early, rapid onset of CT. One 
should expect an earlier onset of CT for meson produc- 
tion than for proton scattering [l4| . as it is much more 
probable to produce a small transverse size in a qq sys- 
tem than in a qqq system. Moreover, the evolution dis- 
tances are easily longer than the nuclear radius , even 
at moderate Q^, which increases the chances for the small 
transverse size object to pass undisturbed through the 
nucleus. These arguments seem to be supported by a 
recent pion-photoproduction experiment at JLab fl3|. 

We report the first measurement of the nucleon num- 
ber, A, and dependence of nuclear transparency for 
the A(e, e'7r+) process. The measurement was performed 
on ^H, ^^C, ^''Al, ®^Cu and ^^''Au nuclei, over a range 
of 1.1 to 4.7 (GeV/c)2. Measurement of both the A and 
dependence of the nuclear transparency is crucial to 
distinguish between CT-like effects and other reaction- 
mechanism based energy dependence of the transparency. 
In this measurement the coherence length for pion pro- 
duction (distance over which the virtual photon fluctu- 
ates into a qq pair) was smaller than the nucleon radius 
and was essentially constant, ranging from 0.2 - 0.5 fm 
over the kinematic range of the experiment. 

The experiment was performed in Hall-C at JLab over 
the summer and fall of 2004. A continuous wave electron 
beam with energies between 4 and 5.8 GeV and currents 



TABLE I: The central kinematics of the experiment. 
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between 10 and 80 /iA was incident on solid foil targets 
of i^C, ^''Al, '^^Cu and ^^''Au, and 4 cm long liquid hy- 
drogen and deuterium targets. The ^''Al target was used 
to mimic the cell walls of the liquid target, and enabled 
the events from the walls to be subtracted from the hy- 
drogen and deuterium data. The scattered electrons and 
the electroproduced pions (tt"*") were detected in coinci- 
dence using the Short Orbit Spectrometer (SOS) and the 
High Momentum Spectrometer (HMS) respectively. The 
kinematic settings of the measurements are shown in Ta- 
ble m The center of mass energy of the hadron system, 
W, for all these settings was above 2.14 GeV in order to 
avoid the resonance region. In addition to the kinematics 
shown in Table U data were collected using the hydrogen 
target over a range of ± 5° around the momentum trans- 
fer direction, 9q = Ott, for each setting (except when 
limited by the hardware constraint 0HMS > 10.6°). These 
data helped develop a model of the elementary pion elec- 
troproduction cross section covering a greater range in 
9q. For two settings (2.15 and 3.91 (GeV/c)^) data 
were also collected at a larger electron scattering angle, 
in order to perform a Rosenbluth separation as a check of 
the reaction mechanism. The results of the Rosenbluth 
separation will be presented in a forthcoming article. 

The SOS gas Cerenkov counter was used to select 
the scattered electrons with an efficiency of better than 
99.2%. The pions were selected using the HMS aero- 
gel [l3| and gas Cerenkov counters, with better than 
98.8% and 99.1% efficiency, respectively. The HMS aero- 
gel Cerenkov counter was used to select pions only for 
HMS central momenta < 3.2 (GeV/c), because they 
were below the pion threshold in the HMS gas Cerenkov 
counter. At higher momentum settings, the HMS gas 
Cerenkov counter alone was sufficient for selecting pions. 
The spectrometer acceptance was determined with a rel- 
ative uncertainty of 1% between targets using a Monte 
Carlo simulation of the experimental apparatus, as de- 
scribed below. For each run, the e— tt'*" coincidence events 
were corrected for random coincidences. The charge 
weighted coincidence yield was also corrected for blocked 
coincidences (< 0.7%), loss of synchronization between 
detectors (< 1.0%), trigger inefficiency ( < 0.5%), elec- 
tronic dead time ( < 1.0%), computer dead time ( < 25%, 
known to much better than 1%), tracking inefficiency ( 
< 4.0%) and particle absorption in the spectrometer ma- 
terial (5.0% known to better than 0.5%). Events from 



3 




FIG. 1: The missing-mass spectra for ^^C(e,e'7r); the crosses 
are data, the solid hne is the simulation which is normalized 
to the data. The shaded region shows the simulated multi- 
pion background. The vertical line indicates the position of 
the multi-pion cut as defined in the text. 

multiplc-pion production were rejected with a cut on the 
missing mass spectrum for each target. Using a simu- 
lation of the multiple-pion background it was estimated 
that the contamination from such events was < 0.4%. 

The standard Hah C Monte Carlo simulation code 
SIMC was used to simulate the experimental appara- 
tus [3]. The p(e,e'7r+)n cross section needed in the 
model was iterated until there was good agreement be- 
tween the simulation and the experimental data. The 
iteration was performed separately for each of the kine- 
matic settings in Table ID A parametrization of the 
p(e, e'7r+)n cross section from previous data fvi\ was used 
as the starting model. 

The Fermi motion of the nucleons in A > 1 targets was 
simulated by folding the elementary cross section with a 
spectral function for the target. For each target, an ap- 
propriate Independent Particle Shell Model (IPSM) spec- 
tral function was used • The simulation includes sev- 
eral corrections such as pion decay within the spectrom- 
eter, external and internal bremsstrahlung radiation and 
pions punching through the collimators at the spectrom- 
eter entrance. It also included corrections due to vari- 
ous reaction mechanism of the A{e,e'TT~^) process such 
as the Coulomb distortion of the incoming and scattered 
electrons, the bound proton in the target nuclei being 
off-shell and Pauli exclusion of the recoiling neutron. 

The phase space for multiple-pion production within 
the spectrometer acceptance was simulated assuming a 
quasi-free single pion production and a uniform phase 
space distribution of the additional pions. The simu- 
lated multi-pion strength was normalized to the tail of 
the experimental missing mass spectra. The multiple- 
pion simulation was used to determine the location of 



the cut on the experimental missing mass spectrum such 
that the contamination from multiple-pion events was 
less than 0.4%. This allowed the missing mass cut to 
be placed ~ 10 - 50 MeV/c^ above the actual kinematic 
threshold for two-pion production. The simulation was 
able to reproduce the shapes of the measured W, and 
\t\ distributions versus the missing mass reasonably well 
for all targets and settings. Representative missing- 
mass spectra for ^^C(e, eV) are shown together with the 
Plane Wave Impulse Approximation (PWIA) simulation 
in Fig. [T] for all settings. The agreement between the 
missing mass spectra obtained from data and simulation 
improve with increasing . The discrepancy seen at 
= 1.1 (GeV/c)^ can be attributed to the reaction mech- 
anisms missing from the simulations such as final state 
interactions between the knocked-out neutron and the 
residual nucleons (nN-FSI) and short range correlations. 

In order to extract the nuclear transparency from the 
experimental yields, the cross section for the bound pro- 
ton must be corrected for the effects of Fermi motion, 
Pauli blocking, the off-shell properties of the proton and 
the acceptances of the spectrometers. In order to account 
for these effects, the nuclear transparency was formed us- 
ing the experimental charge normalized yield, Y, divided 
by the charge normalized Monte Carlo equivalent yield, 
^MC- For a given target, with nucleon number, A, the 
nuclear transparency was defined as: 

T^{Y/YMc)JiY/YMc)^, (1) 

where the denominator is the ratio of the yields from the 
target. As the Monte Carlo simulation does not in- 
clude final-state interactions between the pion and the 
residual nucleons, the nuclear transparency is a mea- 
sure of these final-state interactions, and the reduction 
of these interactions is a signature of CT. 

Traditional nuclear physics calculations based on the 
Glauber multiple scattering mechanism ^1^ are expected 
to be energy-independent ( because the 7r-nucleon cross 
section is constant for the energies in this experiment). 
To investigate the energy dependence, the extracted nu- 
clear transparency is shown as a function of in Fig. [2l 
The point-to-point (Q^ dependent) systematic uncer- 
tainty is 2.4 - 3.2%, dominated by uncertainty in the 
spectral function (1%) and the iteration procedure (1%). 
There is an additional normalization systematic uncer- 
tainty of 1.1% (not shown in the figure) with pion absorp- 
tion correction (0.5%), and target thickness (1%) being 
the main sources. The dependent model uncertainty 
is 7.6%, 5.7%, 3.5%, 3.8% and 3.8% for = 1.1, 2.1, 
3.0, 3.9 and 4.7 (GeV/c)^ respectively. This uncertainty 
was determined from the change in dependence of the 
transparency when using two different spectral functions 
and two different Fermi distributions in the simulation, 
and the dependent uncertainty from reactions mecha- 
nisms not included in the simulation (estimated by quan- 
tifying the difference in shape of the missing-mass spectra 
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FIG. 2: Nuclear transparency, T, vs. for and ^^C (left, 
top panel), ^^Al (right, top), "^^Cu (left, bottom) and ^^''Au 
(right, bottom). The inner error bars are the statistical un- 
certainties and the outer error bars are the statistical and 
point-to-point systematic uncertainties added in quadrature. 
The dark band in the bottom right panel is the dependent 
model uncertainty, and is the same for all nuclei. The solid 
and dashed lines are Glauber and Glauber plus CT calcula- 
tions, respectively [S^]. Similarly, the dot-dash and dotted 
lines are Glauber and Glauber plus CT calculations, respec- 
tively [21!]. These calculations also include the effect of short 
range correlations (SRC). 

from data and simulation) added in quadrature. The 
dependent model uncertainty is shown as a dark band in 
the bottom right panel of Fig. [51 There is an additional 
7.0% normalization type model uncertainty, independent 
of Q^, not shown in the figure. The observed depen- 
dence of the transparency deviates from the calculations 
without CT of Larson et al. and Cosyn et al. (201, i2l!], and 
are in better agreement with the CT calculations of the 
same authors. Larson et al. use a semi-classical Glauber 
multiple scattering approximation, while Cosyn et al. use 
a relativistic version of Glauber multiple scattering the- 
ory. Both groups incorporate CT using the quantum dif- 
fusion model of Ref. [22| with the same parameters r — 
1 fm/c and = 0.7 GeV^. 

In addition to the dependence, the dependence of 
the nuclear transparency on A, is important in the search 
of CT effects and is examined by fitting the transparency 
as a function of A at fixed to the form T = A°'~^. 
The parameter a is found to be ~ 0.76 in fits to the 
pion- nucleus scattering cross sections [23|, and it is ex- 
pected to be energy independent. An energy dependence 
of the parameter a (which quantifies the A dependence 
of nuclear transparency) is a signal for CT-like effects. 
Our results shown in Fig. 3, indicate that the energy de- 
pendence of the parameter a deviates significantly from 
the conventional nuclear physics expectation. The sys- 
tematic uncertainties shown include contributions from 
the fitting error and the model uncertainties. Our re- 
sults are in reasonable agreement with a extracted from 
the calculations (with CT) of Larson et al. [20] but are 
systematically lower than the calculations (with CT and 
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FIG. 3: The parameter a (from T = A°"^) is shown vs g^. 
The inner error bars are the statistical uncertainty and the 
outer error bars are the quadrature sum of statistical and 
systematic and model uncertainties. The hatched band is the 
value of a extracted from pion- nucleus scattering data [23I ]. 
The solid, dashed, and dotted lines are a obtained from fitting 
the A dependence of the theoretical calculations, Glauber, 
Glauber-HCT [l^, and Glauber+SRC+CT HH respectively. 

short range correlations) of Cosyn et al. [2l[. 

These results seem to confirm the predicted early on- 
set of CT in mesons compared to baryons. Our results, 
together with the previous meson transparency measure- 
ments ^lli, il2|] , suggest a gradual transition to meson pro- 
duction with small inter-quark separation, and the onset 
of reaction mechanisms necessary for QCD-factorization 
at values of a few (GeV/c)^. These results also put 
severe constraints on early models of CT which predict a 
dramatic transition with a threshold-like behavior. 

In summary, we have measured the nuclear trans- 
parency of pions from ~ 1.1 to 4.7 (GeV/c)^ over 
a wide range of ^ (2 - 197). Both the energy dependence 
and the A dependence of the transparency show devia- 
tions from the traditional nuclear physics expectations 
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and are in agreement with CT calculations 
is important to extend these measurements to 
(GeV/c)^, where the largest CT effects are predicted, in 
order to establish the onset of CT effect on a firm footing. 
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